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EFFECTS  OF  PARALLEL-JET  MIXING  ON  DOWNSTREAM  MACH  NUMBER 
AND  STAGNATION  PRESSURE  WITH  APPLICATION  TO  ENGINE 
TESTING  IN  SUPERSONIC  TUNNELS 
By  Harry  Bernstein 


SUMMARY 


A one-dimensional  analysis  of  the  results  of  the  parallel-jet 
mixing  encountered  in  the  testing  of  engines  in  supersonic  wind 
tunnels  is  reported.  Equations  were  derived  for  determining  the 
total  pressure  and  Mach  dumber  behind  the  tunnel  terminal  shock. 

The  method  represents  a simple  procedure  for  determining  these 
quantities  while  a tunnel  is  still  in  the  design  Btage.  A specific 
example  of  the  method  is  included. 


INTRODUCTION 


The  results  of  mixing  two  streams  of  different  temperatures, 
pressures,  and  Mach  numbers  are  of  importance  in  the  wind-tunnel 
testing  of  engines.  In  a tunnel  the  exhaust  gases  of  the  engine  mix 
with  the  tunnel  air.  The  effect  of  such  a mixing  upon  the  total 
pressure  and  Mach  number  of  the  combined  streams  is  of  particular 
interest.  The  evaluation  of  these  effects  as  functions  of  engine - 
performance  parameters  is  of  importance  in  determining  the  tunnel 
power  requirements. 


The  present  investigation,  made  at  the  NACA  Lewis  laboratory, 
is  concerned  with  the  analysis  of  the  results  of  this  mixing  process 
by  means  of  one-dimenBional-flow  equations.  The  analysis  is 
restricted  to  the  case  in  which  the  area  of  the  resultant  stream  is 
equal  to  the  sum  of  the  areas  of  the  original  streams.  As  a result 
.of  this  restriction,  the  solution  is  made  dependent  upon  the 
upstream  conditions  only.  The  existence  of  any  area  changes  would 
make  it  necessary  to  evaluate  any  axial  forces  and  include  them 
in  the  momentum  equation.  Such  forces  would  be  dependent  upon  the 
strengths  and  positions  of  any  shocks  in  the  mixing  region;  these 
in  tura  would  be  dependent  upon  the  downstream  static  pressure. 
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The  solution,  when  area  changes  exist,  Is  therefore  dependent,  for  a 
given  set  of  upstream  conditions,  upon  the  downstream  static  pressure. 

The  solution  as  presented  is  also  applicable  to  problems  associated 
with  the  mixing  in  long  cylindrical  ejectors  and  to  any  other  problems 
to  which  the  restrictions  made  in  the  analysis  apply. 


SYMBOLS 

The  following  symbols  are  used  in  this  report : 

A stream  cross-sectional  area 

Cp  specific  heat  of  air  at  constant  pressure 
H heating  value  of  fuel,  Btu/lb 
K ^3/^1 j defined  in  eq.  (9) 

M Mach  number 

air  flow  through  engine,  lb /sec 
nif  fuel  mass-flow  rate,  lb/sec 

mj.  total  mass-flow  rate  of  two  streams  being  mixed,  lb/sec 
(mass  of  fuel  not  included) 

P total  pressure 

p static  pressure 

T total  temperature,  °R 

t static  temperature,  °R 

X defined  in  eq.  (7) 

T defined  in  eq.  (8) 

a product  of  area  and  static  pressure  ratios  of  two  streams  being 

mixed,  A2p2/A1p1 

7 ratio  of  specific  heats  for  air 


M 
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3 


^c 

e 

K 


combustion  efficiency  of  engine, 


AT  across  combustion  chamber 


maximum  possible  AT 
ratio  of  static  temperatures  in  streams  being  mixed,  t2/t-L 

ratio  of  Mach  numbers  of  streams  being  mixed,  Mg/M^ 


t ratio  of  total  temperatures  in  streams  being  mixed,  T2/T]_ 


M* 


a function  of  Mach  number, 


(l  + 

(1  + ?M2)' 


Subscripts : 

m measured  value  (fig.  6) 

max  maximum  measured  value  (fig.  6) 

0 conditions  in  engine  inlet  stream  tube 

1 conditions  in  tunnel  stream  adjacent  to  engine  exit  station 

2 conditions  in  exhaust  stream  at  engine  exit  station 

3 conditions  after  mixing 


Superscripts : 

' conditions  behind  normal  shock  (corresponding  to  unprimed  conditions 
of  same  subscript  ahead  of  shock) 

" conditions  in  resultant  stream  after  heat  addition  in  mixing 

region  (corresponding  to  quantities  of  same  subscript  before 
heat  addition) 


Mims  EQUATIONS 
Derivation 

The  analysis  is  restricted  to  the  case  of  one -dimensional  flow, 
where  the  area  of  the  streams  after  mixing  is  equal  to  the  sum  of  the 
areas  of  the  two  streams  being  mixed  (see  sketch) . The  momentum 
equation  may  contain  no  body  force  terms  between  the  initial  and 
final  stations  at  which  the  flow  is  uniform. 
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Mixing  region 


2 

The  equations  relating  conditions  "before  and  after  the  mixing 
process  are  (since  + Ag  = Aj) : 

Momentum  equation 

PjA^l  + T^!2)  + + ^22)  = + ^32) 

Energy  equation 

PlAlA/^l  (X  + “l2)  + ^ + 2LiJ:  “27 

= P3A3  *Fs  M3  (X  + 

Continuity  equation 

PlAlM!  P2A2M2  _ P5A5M5 

Vtl  V^2 

The  assumption  of  constant  7 and  Cp  was  made  in  order  to 
avoid  the  necessity  of  a trial-and -error  solution.  The  validity  of 
this  assumption  is  discussed  in  appendix  A. 

The  effects  of  viscosity  in  the  "boundary  layer  along  the  wall  of 
the  channel  (friction  force)  are  neglected  throughout  this  analysis. 
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In  most  tunnels  suitable  for  engine  testing,  this  force  is  small 
relative  to  the  pressure  and  momentum  changes  of  the  main  stream. 
The  effects  of  viscosity  and  turbulence  in  the  mixing  process  are 
irrelevant  if  conditions  at  station  3 are  uniform,  as  assumed. 

Division  of  each  of  equations  (l),  (2),  and  (3)  by  its  first 
term  yields 


1 + aX  = 


P3A3  ( 1 + ) 


PlAx(l  + ^M]2) 

P3A3"V^3  M3( 1 + 7 2 1 M32) 

i + o£y  -/e  = — — 4 

PlAiVtl  Mi^l  + 2-^J,  m!2J 


1 + 


qg 

•^/0  PlA^xV^ 


AgP2  Mg 

where  a,  £ , and  0 are  defined  as  -r , — , and 

A^Px  Mx5 


(4) 


(5) 


(6) 


v— , respectively. 
tl 


The  parameters  X and  Y are  defined  as  follows: 


X = 


1 + yg2Mx2 
1 + tMx2 


(7) 


1 + 2-^  ^2mx2 


1 +2L^J, 

2 


“l‘ 


(8) 


Equations  (4),  (5),  and  (6)  may  be  combined  to  eliminate  the 
unknowns  P3  and  tj,  leaving  M3  as  the  only  remaining  unknown. 

Multiplying  equations  (5)  and  (6),  dividing  by  the  square  of 


ana  aiviaing  Dy 

M2  ( 1 + 2 1 M2j 

equation  (4),  and  setting  — L = 

(1  + tM2)2 

(1  + *^(l  + S) 


= <p  yields 


-JeJ  *5  _ 

(1  + oX)2  (Pi 


(9) 


which  shall  be  designated  the  mixing  equation. 
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The  known  conditions  at  stations  1 and  2 and  equation  (9)  are  used 
to  find  the  value  of  f 3.  The  function  cp  is  plotted  against  Mach 

number  in  figure  1.  The  curves  show  that  for  values  of  <p  greater  than 
approximately  0.105,  M is  double-valued  in  <p.  The  two  values  of  M 
are  the  supersonic  and  subsonic  Mach  numbers  on  either  side  of  a normal 
shock.  This  may  be  shown  analytically  from  a solution  of  the  equations 
governing  changes  across  a normal  shock.  As  no  restrictions  were  made 
in  the  derivation,  it  is  inferred  that  (for  the  subsonic  solution) 
the  final  conditions  are  independent  of  the  order  of  occurrence  of  the 
processes  of  mixing  and  shock.  In  addition,  the  existence  of  only 
two  solutions  for  the  Mach  number  after  mixing  implies  independence  of 
the  type  of  shock  system  in  the  mixing  region. 

To  determine  which  value  of  M3  is  physically  correct,  three  cases 
can  be  distinguished: 

(a)  The  mixing  streams  are  both  subsonic.  In  this  case  the  subsonic 
solution  M3 1 is  always  the  correct  one,  since  the  supersonic  solution 

will  result  in  a net  decrease  in  entropy. 

(b)  The  mixing  streams  are  both  supersonic.  In  the  mixing  of  two 
supersonic  streams  both  solutions  are  physically  possible.  The  super- 
sonic solution  will  occur  when  the  conditions  require  that  the  back 
pressure  (static  pressure  after  mixing)  be  low.  If  this  back  pressure 
is  high,  shocks  will  exist  in  the  mixing  region  and  the  resultant 
stream  will  be  subsonic.  Choice  of  the  solution  is  therefore  dictated 
by  the  particular  problem  to  which  the  results  are  applied. 

(c)  One  stream  is  supersonic  while  the  other  is  subsonic.  In  this 
case  the  subsonic  solution  is  always  possible.  The  supersonic  solution 
is  valid  in  cases  where  the  subsonic  jet  is  much  smaller  than  the 
supersonic  stream.  The  existence  of  such  a supersonic  solution  is 
determined  by  an  investigation  of  the  entropy  change  in  the  process. 

If  both  solutions  are  possible,  the  choice  is  again  dependent  upon 
the  particular  problem  under  consideration. 

Tn  many  instances  the  problem  of  mixing  two  streams  of  equal  total 
temperature  (t  = l)  will  arise.  Here,  since  0 = t/y,  equation  (9) 
reduces  to 


When  a solution  for  the  Mach  number  after  mixing  has  been  obtained 
(either  from  eq.  (9)  or  eq.  (10)),  substitution  may  be  made  in  equations  (4) 
and  (6)  to  find  p3  and  t3,  respectively. 
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Discussion 

The  complexity  of  the  mixing  equation  makes  it  difficult  to  determine, 
by  inspection,  the  effects  upon  the  resultant  stream  caused  by  variation 
in  the  values  of  each  of  the  nondimensional  parameters.  In  an  effort  to 
illustrate  some  of  these  effects,  sample  calculations  for  an  of 

3.0  were  performed.  The  results  are  plotted  in  figures  2 and  3.  The 
subsonic  Mach  number  after  mixing  M3 t is  plotted  in  figure  2 as  a 

function  of  the  parameter  a for  a few  chosen  combinations  of  £ and 
0.  The  total-pressure  ratio  Pj./P^  is  plotted  in  figure  3 as  a function 

of  the  same  independent  variables.  The  combinations  of  £ and  0 were 
chosen  to  illustrate  separately  the  effects  of  different  Mach  numbers  and 
different  static  temperatures  in  the  mixing  streams. 

As  seen  in  figure  2,  mixing  streams  of  different  temperatures 
(0  / 1)  caused  an  increase  in  the  subsonic  Mach  number  after  mixing 
over  the  reference  value  (K  = 1,  0=1;  normal-shock  solution  for 
M]^  = 3.0).  This  increase  is  greatest  at  a = 1,  for  any  value  of  0. 

As  a approaches  0 or  «,  the  M3,  curves  approach  the  reference  line 

asymptotically  as  would  be  expected,  since  these  values  of  a imply 
that  the  area  of  one  of  the  mixing  streams  is  negligible  with  respect 
to  the  other. 

The  effects  of  mixing  streams  of  different  Mach  number  (£  ^ l) 
are  observed  to  vary  greatly  with  the  value  of  £;  however,  a few 
general  conclusions  about  these  variations  may  be  stated.  For  any  value 
of  £ , the  curve  is  asymptotic  to  the  reference  line  as  a approaches 
zero.  If  a value  of  £ > 1 is  used,  the  resulting  curve  rises  to  a 
peak  at  some  value  of  a less  than  unity.  The  curve  then  drops  below 
the  reference  line  and  approaches  a value  of  M3 . asymptotically  as 

a approaches  «*.  This  value  of  M3,  is  the  subsonic  Mach  number 
behind  a normal  shock  at  a Mach  number  M2  (M2  = £M^) . For  § < 1, 

the  curves  rise  above  the  reference  line  and  reach  a peak  at  some 
value  of  a larger  than  those  values  shown  on  the  figure.  (Differ- 
entiation of  the  mixing  equation  with  respect  to  a shows  that  the 
curve  of  resultant  Mach  number  has  one  peak  for  any  values  of  £ and 
0,  except  the  combination  £ = 1,  0=1.)  As  a approaches  ® these 
curves  are  also  asymptotic  to  some  Mach  number  M31.  This  value  is 

either  the  subsonic  Mach  number  behind  a shock  at  Mach  number  M2 
(for  M2  > 1) , or  M2  itself  (for  M2  < l) . 

In  the  computations  for  the  total-pressure  ratios  P3 , /P^,  some 

assumption  had  to  be  made  concerning  the  variable  a.  While  the  Mach 
number  curves  (fig.  2)  are  quite  general,  the  total-pressure-ratio 
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curves  depend  upon  the  actual  components  of  a,  that  Is,  upon  Pg/Pi  811(1 
Ag/Ai  (see  eq.  (4)).  The  ratio  Pg/Pi  was  therefore  arbitrarily  chosen 
to  be  unity,  and  hence  a,  in  figure  3,  represents  the  area  ratio  Ajj/a^ 

For  4=1  and  different  static  temperatures  (0  ^ 1),  the  pressure- 
ratio  curves  (fig.  3(a))  fall  below  the  reference  line,  become  asymptotic 
to  this  line  as  a approaches  0 or  and  have  a minimum  point  at 

a = 1.  The  decrease  in  total  pressure  is  due  to  the  heat  exchange  between 
the  streams  being  mixed  and  the  net  entropy  increase  associated  with  this 
heat  exchange. 

For  4 £ 1 (unequal  Mach  numbers),  the  curves  (fig.  3(b))  are 
asymptotic  to  the  reference  line  as  a approaches  zero.  As  a increases, 
the  curves  depart  from  the  reference  line,  falling  above  or  below  it  as 
4 is  greater  than  or  less  than  unity,  respectively.  As  a approaches 
infinity,  these  curves  approach  the  values  Pjj./p^  (for  M2  = 4M^  > l) 

or  P2/P-l  (for  M2  < !•)  • 


APPLICATION  TO  ENGINE  TESTING  IN  SUPERSONIC  WIND  TUNNELS 

The  equations  just  derived  may  be  applied  in  the  design  of  a super- 
sonic wind  tunnel  in  which  engines  are  to  be  tested.  In  such  a tunnel 
the  exhaust  gases  of  the  engine  mix  with  the  tunnel  air;  this  mixing 
affects  the  values  of  the  resultant  Mach  number  and  total  pressure.  The 
evaluation  of  the  resultant  stream  properties  as  functions  of  engine- 
performance  parameters  is  of  importance  in  determining  the  power  and 
pres sure -ratio  requirements  of  the  tunnel.  The  arrangement  is 
schematically  illustrated  in  figure  4. 

At  the  engine  exit  station,  mean  values  of  the  flow  properties  of 
the  tunnel  stream  (station  l)  are  required  for  the  one -dimensional 
analysis.  If  A2/Aq  1b  close  to  unity  and  the  engine  is  operating  at  a 

mass-flow  ratio  of  unity,  these  mean  conditions  may  be  taken  as  free- 
stream  conditions.  More  accurate  mean  values  may  be  obtained,  if 
necessary,  by  constructing  the  flow  field  past  the  outer  surface  of  the 
engine . 

Flow  properties  in  the  exhaust  jet  at  the  engine  exit  (station  2) 
may  be  evaluated  as  functions  of  engine -performance  parameters  (P2/Pq, 

M2,  t).  Again,  the  assumption  of  an  engine  operating  at  a mass-flow 

ratio  of  unity  is  made.  If  the  mass  of  the  fuel  is  neglected,  the 
condition  of  constant  mass  flow  yields  the  following  relation  between 
engine  pressure  recovery,  exhaust- jet  Mach  number,  and  total-temperature 
ratio: 
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G + ^ir1  “a2) 

_ ^2  ^2 

G + Mo2 

f 

L m2 

P0  A0 

Mo 

_ 

1 


(11) 


If  the  geometry  of  the  engine  to  be  tested  (M2  is  a function  of  exhaust- 

nozzle  geometry  only)  and  t are  known,  equation  (ll)  provides  a means 
of  calculating  the  total  pressure  in  the  exhaust  jet. 


Once  the  conditions  in  the  exhaust  jet  and  the  adjacent  tunnel 
stream  have  been  determined,  use  may  be  made  of  the  mixing  equation 
(eq.  (9))  to  evaluate  tunnel  Mach  number  and  total  pressure  after 
mixing.  For  this  application,  station  3 must  be  assumed  to  be  a suitable 
distance  downstream  of  the  tunnel  terminal  shock,  and  hence,  the  subsonic 
solution  to  the  mixing  equation  is  the  proper  solution.  The  reasons  for 
this  assumption  are  clearly  illustrated  in  figures  5 and  6.  Figure  5 
shows  a jet  exhausting  into  a supersonic  stream.  The  jet  is  observed 
to  expand  slightly  to  satisfy  ambient  static -pressure  conditions,  but 
little  mixing  is  observed.  This  is  further  evidenced  by  the  temperature 
profiles  presented  in  figure  6.  These  profiles,  which  are  affected  by 
the  amount  of  mixing,  show  that  the  major  portion  of  the  mixing  occurs 
in  a region  downstream  of  the  tunnel  terminal  shock,  and  hence  the 
requirement  that  station  3 be  a suitable  distance  downstream  of  the 
terminal  shock. 

Additional  changes  in  Mach  number  and  total  pressure  may  be  caused 
by  heat  addition  in  the  mixing  region.  This  is  usually  the  case  in 
engine  testing,  for  excess  fuel  Is  carried  out  of  the  engine  in  the 
exhaust  jet.  This  fuel  bums  in  the  mixing  region  just  downstream  of 
the  tunnel  terminal  shock.  Figure  7 illustrates  this  effect.  Inasmuch 
as  little  fuel  is  usually  carried  out  with  the  exhaust  gas,  as  compared 
with  the  total  mass-flow  rate  mt,  the  magnitude  of  any  changes  due 

to  this  burning  may  be  small,  and  in  many  cases  can  be  neglected. 

In  appendix  B,  relations  are  derived  which  may  be  used  in  evaluating 
the  changes  due  to  thiB  heat  addition,  if  a higher  degree  of  accuracy 
is  desired. 

Throughout  this  analysis  the  tunnel  terminal  shock  was  assumed  to 
be  located  in  the  test  section  downstream  of  the  model.  This  represents 
peak  efficiency  operation  of  a tunnel  having  no  second  throat.  While 
this  is  possible,  most  wind  tunnels  are  not  operated  at  peak  efficiency 
but  are  operated  with  the  terminal  shock  positioned  a short  distance 
downstream  of  the  start  of  the  tunnel  diffuser.  If  this  is  the  case, 
the  results  of  this  peak  efficiency  analysis  must  be  corrected  for  the 
additional  diffuser  losses.  In  general,  however,  the  trends  indicated 
by  this  analysis  will  be  unaffected.  This  constant-area  analysis  is 
not  applicable  for  tunnels  with  second  throats. 
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EXAMPLE 

Consideration  has  been  given  to  the  problem  of  testing  an  engine 
having  an  8 -inch  outlet  diameter  in  an  18-  by  18 -inch  tunnel  operating 
at  a Mach  number  of  3.1.  To  be  found  are  the  effects  upon  tunnel 
operating  conditions  caused  by  variations  in  t and  the  exit- nozzle 
throat  area  of  the  engine  being  tested. 

In  this  solution,  the  engine-air  mass-flov  ratio  was  assumed  to  be 
unity.  Equation  (11)  was  used  to  compute  values  of  engine  pressure 
recovery  P2/Pq  ^or  various  values  of  exit  Mach  number  M2  and  v . 

In  order  to  simplify  the  problem,  it  was  assumed  that  A2/Aq  = 1 and 

that  the  engine  was  a perfect  cylinder,  making  the  conditions  at 
stations  0 and  1 equal  (see  fig.  4). 

With  conditions  at  stations  1 and  2 known,  the  conditions  in  the 
tunnel  after  mixing  were  obtained  by  use  of  equations  (9)  and  (4).  The 
results  of  these  computations  are  plotted  in  figures  8 and  9. 

Figure  8 shows  tunnel  Mach  number  after  mixing  downstream  of  the 
normal  shock  in  the  test  section  M3  ( as  a function  of  engine  pressure 

recovery  for  various  total  temperature  ratios.  Lines  of  constant  engine- 
exit  Mach  number  are  shown.  Figure  9 shows  resulting  tunnel  pressure 
recovery  P3 > /P^  as  a function  of  the  same  independent  variables. 

Increases  in  total-temperature  ratio,  while  holding  either  M2  or 
P2/Pq  constant,  result  in  higher  values  of  M3 1 and  P3t/P-^.  An 

increase  in  the  engine  pressure  recovery,  while  holding  the  total- 
temperature  ratio  constant,  increases  the  tunnel  pressure  recovery  and 
decreases  the  value  of  M3 . . 

The  increase  in  the  tunnel  pressure  recoveries  with  increasing 
total -temperature  ratios  (greater  heat  addition)  is  explainable  as 
follows.  From  the  theory  of  one-dimensional  gas  flow  (ref.  l),  these 
facts  may  be  stated  about  conditions  at  station  3':  (l)  At  a constant 

value  of  total  temperature,  decreases  in  total  momentum  (pA(l  + }M2)) 
at  this  station  will  result  in  increased  values  of  M3i  and  decreased 

tunnel  pressure  recovery  P3I/P1,  and  (2)  If  total  momentum  is  held 
constant,  increases  in  energy  will  result  in  greater  values  of  M3 1 
and  smaller  tunnel  recoveries. 

In  figures  8 and  9,  lines  of  constant  t are  lines  of  constant 
energy  (constant  T3,).  Lines  of  constant  total  momentum 

(aX  *=  constant)  have  been  drawn  on  each  figure.  Larger  values  of 
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aX  are  synonomous  with  greater  total  momentum  at  station  3’  (eq..  (4)), 
as  station  1 represents  tunnel  free-stream  conditions  which  are  constant. 
If  changes  are  observed  along  these  energy  and  momentun  lines,  the 
results  are  seen  to  he  in  agreement  with  the  one -dimensional  gas-flow 
theory.  Variations  in  Pji/P^  with  changes  in  fuel-air  ratios  are 

therefore  an  effect  of  combined  momentum  and  energy  changes. 


CONCLUDING  REMARKS 

A one-dimensional  analysis  of  the  results  of  the  parallel- jet 
mixing  encountered  in  the  testing  of  engines  in  supersonic  wind  tunnels 
has  been  reported.  This  type  of  analysis  presents  a reasonable  approach 
to  obtaining  approximate  figures  for  the  tunnel  operating  conditions 
while  the  tunnel  is  still  in  the  design  stage.  These  figures  would  be 
based  upon  the  known  tunnel  geometry  and  inlet  conditions  and  estimations 
of  the  model  geometry  and  values  of  the  engine -performance  parameters. 
Additional  equations  are  presented  for  evaluation  of  changes  due  to  the 
burning  of  excess  fuel  downstream  of  the  engine -exhaust  station. 

In  the  evaluation  of  the  properties  of  the  engine -exhaust  jet  and 
the  resultant  (mixed)  stream,  it  has  been  demonstrated  that  the 
assumptions  of  constant  Cp  and  y Introduce  no  significant  errors  in 

the  final  results. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  January  5,  1953 
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APPENDIX  A 

DISCUSSION  OF  ERROR  INVOLVED  IN  ASSUMPTION  OF  CONSTANT  y AND  Cp 

In  view  of  the  values  of  temperature  associated  with  the  testing  of 
engines  in  supersonic  wind  tunnels  and  the  variations  of  Cp  and  y at 

these  temperatures,  the  assumption  of  constant  Cp  and  y,  as  made  in 

the  mixing  solution  and  in  the  determination  of  exhaust- jet  properties, 
appears  justified. 

To  illustrate,  it  is  assumed  that  the  total  temperature  of  the 
exhaust  Jet  is  limited  to  a maximum  of  3000°  R,  a value  based  approx- 
imately on  the  highest  temperatures  which  present-day  materials  can 
withstand.  It  is  also  assumed  that  the  exhaust -jet  Mach  number  will  be 
greater  than  or  equal  to  1.6,  a reasonable  value  for  an  engine  operating 
in  a supersonic  stream  of  moderate  Mach  number.  These  factors  place 
an  upper  limit  of  approximately  2000°  R on  the  static  temperature  of  the 
exhaust,  stream,  for  an  extreme  case.  It  is  to  be  realized  that  as  the 
exhaust-jet  Mach  number  increases  or  total  temperature  decreases  or  both, 
the  value  of  its  static  temperature  decreases. 

In  the  usual  supersonic  wind  tunnel,  air  is  expanded  to  a tempera- 
ture of  the  order  of  200°  S,  although  this  figure  may  be  decreased  con- 
siderably for  hypersonic  tunnels.  Hence,  the  range  of  static  temperatures 
of  importance  in  the  mixing  problem  is  from  200°  to  2000°  R.  The  mixed 
stream,  prior  to  the  tunnel  terminal  shock,  would  have  a static  tempera- 
ture greater  than,  but  much  closer  to,  the  200°  R figure,  if  supersonic 
mixing  existed.  This  is  a result  of  the  small  mass  flow  through  the 
engine  (hot  air)  as  compared  with  that  through  the  tunnel  adjacent  to 
the  engine.  The  tunnel  terminal  shock  may  cause  the  static  temperature 
to  increase  by  a factor  of  2 or  3,  but  a value  of  2000°  R is  still  felt 
to  be  an  upper  limit  of  static  temperature  in  a very  extreme  case. 

For  air,  the  variations  in  Cp  and  y are  (ref.  2): 


°R 

V 

Btu 

lb-°R 

7 

200 

0.2395 

1.400 

500 

.2400 

1.400 

1000 

.2488 

1.380 

1500 

.2644 

1.349 

2000 

.2776 

1.328 

2418 
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For  2000°  R the  variations  in  Cp  and  y from  the  values  of  200°  R 
are  15.8  and  5.2  percent,  respectively,  although  these  variations  decrease 
rapidly  for  temperatures  less  than  2000°  R.  In  the  range  of  temperatures 
from  200°  to  540°  R,  hoth  Cp  and  y are  constant. 

To  show  the  effects  upon  the  accuracy  of  the  results  of  the  mixing 
problem  when  y and  Cp,  are  assumed  constant,  the  following  problem 
was  considered : 

Tunnel  size,  18  in.  by  18  in. 

Model  size,  8 in.  diam. 

M0  = 3.1 

M2  = 1.76 
T0  = 550°  R 
T2  = 2370°  R 
t)c  = 100  percent 

The  exhaust- jet  properties  were  found  from  these  data,  and  then  the 
mixing  equations  were  applied,  both  with  Cp  and  y assumed  constant 

and  with  values  of  Cp  and  y dependent  upon  the  temperature.  The 

latter  method  involved  a trial-and-error  solution,  the  details  of 
which  are  not  presented  here.  A tabulation  of  results  for  both  cases, 
along  with  percentage  variations,  follows: 


cp,2 

Variation 
in  Cp, 

percent 

?2 

Variation 
in  7, 
percent 

Variation 
in  M3 . , 

percent 

P3'/Pl 

Variation 

in 

P3'/pi> 

percent 

Constant  Cp 
and  7 

0.2400 

9.6 

1.400 

3.4 

0.580 

2.6 

0.3326 

1.1 

Variable  Cp 
and  7 

0.263 

1.352 

0.565 

0.3288 

large  variations  in  Cp  and  y are  seen  to  result  in  very  small  varia- 
tions in  the  final  results. 
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APPENDIX  B 

HEAT  ADDITION  IN  MIXING  REGION 


The  equation  governing  the  flow  of  a fluid  in  a constant-area  channel 
with  heat  addition  can  he  written 

Continuity  equation 


pM  ff.  + M2"  + rZI  m"2 

Momentum  equation 

p(l  + ?M2)  = p"(l  + 7M"2) 


(Bl) 


(B2) 


Division  of  the  continuity  equation  by  the  momentum  equation  yields 


M fj  1 + 2Lf^-  M2 

r^r 

M"  ljl  + 7-~  -1-  M"2 

(1  + tm2) 

j T ~ 

(1  + rM"2) 

(B3) 


Equation  (B3),  when  squared,  becomes 


(B4) 


Application  of  this  relation  to  the  solution  of  the  mixing  problem 
yields 


(p3 


<p3,t 


(B5) 


Therefore  it  is  seen  that 
the  equation 


q>3it  can  be  directly  obtained  from 


*1  fcy 


(B6) 


The  ratio  T3"/t3  is  given  in  terms  of  fuel  heating  value  and 
engine  combustion  efficiency  as 


Su 

i**Jlw  ‘ *•• 
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M 

X> 

P- 


& 

| 


I 

| 

I 

! 


T3” 

t3 


(1 


= 1 + 


- '•{?> 


CpT3 


(B7) 


where 


^f  _ mf  ®a  ^1  _ JJJf  AxgV  1 


“t  “a  ml  mt  ma  V 1 + StL 

Vey 


(B8) 


and 


'1  + 


t3  = T1 


ogT 


1 + 


gg 

Ve- 


Combining  equations  (B7),  (B8),  and  (B9)  yields 


(B9) 


(BIO) 


The  combustion  efficiency  tj  is  given  approximately  as 

c 


= 


_ (T  ~ 3-)cpTi 


mj. 

mD 


H 


(Bll) 


When  <P3»  has  been  evaluated,  M3"  may  be  obtained  from 

figure  1.  Choice  of  the  subsonic  or  supersonic  solution  is  determined 
by  the  same  factors  discussed  in  relation  to  equation  (9)  or  supersonic 
wind  tunnels. 

The  conservation  of  momentum  and  mass  flow  during  the  heat  addition 
process  allows  the  use  of  equations  (4)  and  (6)  for  determining  p3n 

and  t3>.,  respectively,  from  the  value  of  M3». 

A comparison  of  equations  (9)  and  (B6)  illustrates  the  effect  of 
heat  addition  in  the  mixing  region.  Since  T3"/T3  is  always  greater 

than  unity  for  such  a heat  addition,  then  <p3.t>  (p3,  and  hence 
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M3t'>M3i.  Addition  of  heat  in  the  mixing  region  therefore  causes  an 

increase  in  the  subsonic  Mach  number  at  station  3.  A further  effect, 
as  indicated  in  reference  1,  is  to  cause  an  increase  in  the  entropy  of 
the  stream,  with  a subsequent  decrease  in  total  pressure  at  station  3. 
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Figure  2.  - Effect  of  variation  in  parameters  of  mixing  equation  on  Mach  number  after  mixing  (M. 


Figure  3.  - Effect  of  variation  of  parameters  in  mixing  equation  on  total-pressure  ratio  of  mixing  process 
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Figure  5 . - Photograph 
supersonic  stream  of 


of  5/8- inch-diameter  jet  exhausting 
Mach  number  2.07. 
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Figure  6.  - Temperature  profiles  for  jet  shown  in  figure  5. 
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(a)  Normal  mixture  (m  /m  = 0.033). 

I Q 
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v — Upstream  edge  of  terminal  shock 
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Figure  8.  - Tunnel  Mach  number  after  mixing  as  a function  of  engine  pressure  recovery  and 
total  temperature  ratio  for  constant  engine  mass  flow.  M-^  = 3.1;  Ag/A-^  = 0.1857;  T = 1 
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Figure  9.  - Tunnel  pressure  recovery  as  a function  of  engine  pressure  recovery  and  total 
temperature  ratio  for  constant  engine  mass  flow.  = 3.1;  ^ = 0.1837;  f = 1.4. 


